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REMARKS 

Claims 1-9 were pending in the application and have been canceled. New claims 10-22 
have been added. Accordingly, claims 10 - 22 are currently pending in the application and are 
presented for reconsideration in view of the following remarks. 

In the Office Action, the Examiner objected to the Abstract because of its undue length and 
grammatical errors; objected to the Specification for not containing the proper layout and section 
headings and because of grammatical errors; objected to claims 1-9 because they are generally 
narrative and indefinite; objected to claims 1-9 because of informalities; and allowed claims 1 -9 if 
rewritten to overcome the claim objections. 

By this Amendment, claims 1 - 9 are canceled and replaced with claims 10 - 19; and new 
claims 20 - 22 have been added. The amendments to claims 10-18 were made to overcome the 
Examiner's objections. Support for new claim 1 9 can be found for example in Figure 2; new claim 

20 can be found for example on page 24 of the clean copy of the substitute specification; new claim 

21 can be found for example in Figure 3; and new claim 22 can be found for example on page 21 of 
the clean copy of the substitute specification. 

It is further respectfully submitted that the within amendments introduce no new matter 
within the meaning of 35 U.S.C. §132. 

Objection to the Abstract 

The Examiner objected to the Abstract because of its undue length and grammatical errors. 

In response, Applicant has amended the Abstract to contain 1 50 words or less and to correct 
grammatical errors. Therefore, Applicant respectfully requests that the objection to the Abstract be 
withdrawn. 

10 



Application No. 10/649,941 

Art Unit: 3663 

Attorney Docket No. 25585 

Reply to Office Action mailed February 25, 2005 

Objection to the Specification 

The Examiner objected to the Specification for not containing the proper layout and section 
headings and because of grammatical errors. 

In response, Applicant has amended the Specification as provided in 37 C.F.R. 
1.77(b) regarding the layout and section headings and to correct grammatical errors. Therefore, 
Applicant respectfully requests that the objection to the Specification be withdrawn. 

Claim Objections 

The Examiner objected to claims 1-9 because they are generally narrative and indefinite; 
objected to claims 1-9 because of informalities; and allowed claims 1 -9 if rewritten to overcome 
the claim objections. 

In response, Applicant wishes to thank the Examiner for the indication of allowable subject 
matter. Applicant has canceled claims 1-9 and replaced them with new claims 10-18. Further new 
claims 19-22 have been added. Applicant respectfully submits that the claims now comply with 
U.S. practice and 35 U.S.C. 112. Therefore, Applicant respectfully requests that the objection to the 
claims be withdrawn and that the Examiner pass the application to issue. 
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Conclusion 



In light of the foregoing, Applicant submits that the application is now in condition for 
allowance. If the Examiner believes the application is not in condition for allowance, Applicant 
respectfully requests that the Examiner contact the undersigned attorney if it is believed that such 
contact will expedite the prosecution of the application. Favorable action with an early allowance of 
the claims is earnestly solicited. 



Respectfully submitted, 



NATH & ASSOCIATES PLLC 



April 25,2005 




NATH & ASSOCIATES PLLC 
1030 15 th Street, N.W. 
6 th Floor 

Washington, D.C. 20005 
Tel: (202)775-8383 
Fax: (202)775-8396 



Reg. No. 26,965 
Gregory B. Kang 
Reg. No. 45,273 
Teresa M. Arroyo 
Reg. No. 50,015 



Customer No. 20529 
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p"g>s. A method of wave diagnostics of the oil - and - gas deposit 
•§>\ BACKGROUND OF THE INVENTION 

& The given invention refers to the area of the exploration geophysical methods ef 
texploration , in particular, to the seismic survey of the oil-and-gas deposits on the basis of the 
complex of waves of the different types. The invention is intended for diagnostics of the reservoir 
and fluid-saturated properties of the rocks, situated in the vicinity of the borehole^ and Further 
diagnostics include fe f the reception of the quantitative estimates of the significant oil-field 
parameters for the identification of the image of the oil-and-gas deposit, and fef the subsequent 
effective baring of the productive layers in the time of exploration and development drilling of the 
boreholes as well as for the usage of the obtained data for the calculation of the oil and gas reserves. 

The method is especially effective for diagnostics of the fissured reservoir (limestones, 
sandstones) because it allows to m e asur e the measurement of such parameters of the fissured rock 
as the average length and linear density of fractures, fissured porosity and permeability, fissured 
residual water saturation and fluid saturation, while the majority of the last cannot be determined 
with the use of the standard complex of the well logging methods. Besides, the method can be used 
for the investigation of the lithological and petrophysical parameters of the rocks in situ. 

The method of the vibroseism exploration of the oil-and-gas deposits (Ru, Patent N° 
2045079 GO IV 1/00, 1995) is known, according to which they generate the seismic oscillations 
with the frequency of £2 - 20) Hz, register the seismic responses of the earth before and after the 
excitation of oscillations on three components not less than two geophones simultaneously, and on 
sharp magnification of the extreme of the received amplitude-frequency characteristics of the 
seismic response at the frequencies of (2 - 6) Hz make the decision about the presence of the oil- 
and-gas deposit. The method allows to det e ct the detection of the presence of the oil-and-gas 
deposit in some cases, though its efficiency is not high as far as the amplitude extreme is the 
parameter depending on the different factors. In particular, the amplitude values essentially depend 
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on the conditions of the excitation and registration, the absence of data about which can give the 
false information about the presence of the oil-and-gas deposit. 

Therefore, for the identification of installation on such the parameter as the amplitude 
extreme it is necessary to install the level of amplitude discrimination (or threshold), and for the 
correct choice of the threshold magnitude the additional information on the statistical properties of 
useful signal and noise is r e qui e r e d required together with the distribution of probabilities of the 
amplitudes excess of the useful signal and noise of various levels. The method does not envisage 
the control over the change of the conditions of the wave generation, and also demands the major 
expenditures of the time at the usage of additional devices, as existing vibrators cannot generate the 
seismic waves in the infrasonic range with the required frequency of the oscillations excitations. 

The method of the investigation of hydrocarbon deposits on the basis of the effects of their 
nonideal elasticity is known (Rapoport M.B., Ryzhkov V.I. 1992-2000. Scientific Reports at 
International Geophysical Conferences of SEG, IGRC, EAGE etc). In the above-mentioned method 
the experimental fact is used, that the hydrocarbon-saturated porous medium has tall absorption 
properties. It is supposed, that the effects of seismic nonelasticity of the rocks are shown by means 
of frequency - dependent absorption and dispersion of the phase velocities of the seismic waves, the 
values of which are calculated according to the data of the surface and borehole seismic survey. On 
the basis of method offered by the authors, the investigations of seismic absorbtion according to the 
surface-seismic survey data (of 2D/3D) were implemented, moreover for the analysis of the 
absorption parameters the nonstacking seismic data and traceability of the seismic reflections on the 
time cross sections were used. The received values of the absorption parameters are used as the 
indicators of the hydrocarbon presence in the surface seismic survey. The vertical seismic profiling 
(VSP) data in the productive boreholes are regarded as parametric at the analysis of seismic 
nonelasticity and are applied for the separation of the productive thicknesses skipped in the time of 
well logging, while the borehole observations of VSP, the abnormality absorption parameters and 
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velocity dispersions received by the surface seismic survey, are immediately connected with the oil- 
and-gas deposits. For the definition of the threshold level, which corresponds to the outline of the 
deposit, the values of the indicator parameters are calibrated permanently over the results of the 
tests of the productive boreholes.The usage of the given method allows to define the most 
perspective in the oil-and-gas-bearing relations of the investigated area, and at the presence of the 
boreholes also allows to delineate the productive sections. 

However, the above-mentioned method uses as calculated of the surface seismic survey data 
and includes the condition of the mandatory adapting of the investigation results with the well 
logging data. Because of lapses in the time of parametric measurement and in the presence ef only 
the effective attenuations estimates of e ff e ctiv e att e nuations of the amplitudes method can give the 
erroneous prognosis concerning the discovery of layers promising in the oil-and-gas-bearing 
relation. 

Most close according to the technical essence to the declared invention is "The-A method of 
the determination of the rheological properties of solid-liquid mediums by means of the multiwave 
borehole seismic survey" (UA, Patent Ns 718 CI G Ol V 1/40, G 01 V 1/00. 15.12.93. Bull, tfo 2), 
implementing the principles of borehole seismic survey (prototype). The essence of the invention 
consists in the following. In the vicinity of the borehole the flat -plane longitudinal wave is formed 
by means of the e xcitation of seismic oscillations excitation by the near-surface sources. By means 
of the three-component geophones of the displacements (or stresses) the dynamical characteristics 
of the direct seismic waves in the investigated borehole are recorded, at the same time one of the 
geophones is oriented along the normal line to the plane wavefront of the longitudinal wave. Then, 
according to the received kinematic characteristics of the longitudinal and shear (or converted) 
waves the numeral processing of the dynamical characteristics of the direct longitudinal wave is 
fulfilled. For this purpose they execute the separation of the monotype t e mporal time signals for the 
registered downhole borehole components of displacements (or stresses) by pointwise deduction out 
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of the axial component (z- component) and the tangential component (x,y- components) of the 

displacements (or stresses) with the obtaining of the significant parametric estimates of the separate 
signals in the future. 

After that according to the known analytical relationships the kinematic and dynamical 
parameters are transformed in respect to the bedded values with the further determination of the 
rheological parameters of the rocks, as well as the absorption parameters of the longitudinal wave 
for each observation point along the depth of the borehole. The known method of the borehole 
seismic survey allows over the dynamical characteristics of the direct (downgoing) seismic waves 
to define the rheological and absorption parameters in porous - fissured layers, which are opened by 
means of the exploration borehole, and in view of the essential differentiation of the oil-and-gas 
deposits over rheological properties the reliable information about the physico-mechanical 
parameters of the rocks in the vicinity of the borehole is received. 

However^ the above-mentioned method ignores the conditions changes at the multiple 
excitation of the oscillations by the near-surface sources and the influence of the specified changes 
on the dynamical parameters of longitudinal wave. The method also does not envisage the spacial 
statistical processing of the amplitude values of the monotype time signals for the receiving of the 
coefficients of attenuation and energy dispersion of longitudinal wave. The ultimate aim of the 
known method does not allow to receive the bedded values of the oil-field parameters for the exact 
identification of the oil-and-gas deposits images . 

The kernel of the problem which advantaged the creation of the present invention, consists 
in the following. The features of the existing technology of the deep borehole drilling with the 
considerable differentiation of bedded pressures in the rocks along depth of the borehole is 
demanded in order the above-mentioned drilling mud heavers will be applied under the specified 
conditions for the prevention of emergencies, therefore, there is the peening of the part of 
perspective in the oil-and-gas-bearing relation layers. By the appearance of the drilling mud in last 
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considerable invaded zones qualitative diagnostics of the productive layers by means of the 
standard methods of the well logging becomes ineffective, what in and the time of the deposit 
development leads to the skip of the average-and-low-permeable oil-and-gas-bearing layers and 
decreases the potential productivity of the boreholes in (1^2 - UJ>) times. 

The task, set at making the given invention, is to receive the method of wave diagnostics of 
the oil-and-gas deposit which by means of the seismic insonification of the vicinity of the borehole 
will allow over the recording of the displacement vector data in the investigated borehole by means 
of the three-component seismic profiling to receive the oil-field information on the values of the 
absorption coefficients, characterizing the type of the fluid saturation; the coefficients of Assuring, 
granular and fissured permeability, residual water saturation and fluid saturation of the rocks, 
including the layers with the considerable invaded zone of the drilling mud. 

The existing problem of transition from the observable values of the parameters to their 
bedded values can be solved by the data special processing, introducing the corrections to the 
dynamical parameters of the downhol e borehole and check data, taking into account the changes of 
the conditions of the wave excitation and the filtrational effect of the covering thickness of the rocks 
for each observation point along depth of the investigated borehole. 

Th e last Last allows first to calculate the dynamical parameters of the impulse responses of 
the layers, and then to transform their corresponding petrophysical and oil-field parameters and 
simultaneously provides the reception of the reliable information about the effective bed-thickness of 
the productive layers and about th e location sites of the gas-oil-water contacts, and also enables to 
estimate the quality of the cap rocks and to establish the presence of the tectonic dislocations of the 
investigated borehole of crossing. 

BRIEF SUMMARY OF THE INVENTION 

The problem is solved in the following way. A method of wave diagnostics of 
the oil-and-gas deposit which includes the excitation of the seismic waves by the near- 
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surface sources of oscillations, the recording of the displacement vectors for the direct longitudinal 
and shear waves for the each observation point by means of the three-component seismic profiling 
along the depth of the investigated borehole and the processing of their kinematic and dynamical 
parameters, characterised that one near-surface source of oscillations executes in the 
vicinity of the borehole the longitudinal wave excitation , and their reception is carried out 
simultaneously in the investigated and check boreholes; after th e spac e orientation in space of data 
of three-component obs e rvation data, observations, the monotype t e mporal time signals for the 
compression and shear constituents of the direct longitudinal wave are extracted from out of the 
seismic signals wavelets which were registered on the axial and tangential components of the 
displacement vector fer by the downhole borehole and check devices respectively; further by means 
of the-processing ef-the kinematic parameters of the direct longitudinal and shear waves the 
calculate interval values of the e lastic moduli elasticity modules (longitudinal and shear) are 
calculated ; after the decoding of the monotype temporal time signals quantitative estimates of their 
significant dynamical parameters are made for e ach component and th e obs e rvation point 
according to the downhol e borehole and check data, and respectively for the each observation point 
and constituent of the direct longitudinal wave, moreover the accuracy of the obtained estimates of 
the dynamical parameters is monitored by means of the computer modelling of the seismic signals; 
Then-then, carrying out the correction of the corresponding parameter because of the changes of the 
wave excitation conditions during multiple excitation of oscillations and the filtration of the seismic 
signals in the covering thickness of the rocks ^th e y d e fin e defines the dynamical parameters of the 
impulse responses, the-values of the coefficients of dynamical viscosity (volumetric and shear) and 
the-absorption (for the-compression and shear compon e nts constituents) of the layers, when the 
calibrated values of the absorption coefficients for the compression constituent are used a s th e 
indicators of the hydrocarbon pr e sence in the capacity of the indicator for the presence estimate and the 
fluid-type in th e layers, identification, filling the pore and fissure volume in the rocks, and also &e-values 
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of the bedded coefficients of the amplitud e effective attenuation, attenuation and energy dispersion 
of the longitudinal wave are calculated; then, using th e functional r e lationships, which ar e taking 
into account th e th e rmodynamical conditions of th e rock b e dding, the obtained petrophysical data 
are converted into the bedded values of the parameters of the reservoir and fluid-saturated 
properties of the layers with the aim of the receipt of the necessary totality of the-significant oil- 
field parameters for the identification of the oil-and-gas deposit image. 

At the same time the d e pth of th e location of th e e xplosive charg e (or air gun) for th e 
e xcitation of oscillations in th e n e ar surfac e check borehole is d e fin e d out of th e condition of th e 
separation of th e is located on the trajectory of propagation direct longitudinal (or shear wave) 
between investigated borehole and the point of the oscillations excitation at the distance (20 - 50) m 
from the last. Location depth of the explosive charge (or air gun) used for the oscillations excitation 
in the near-surface borehole is defined out of the separation condition of the direct longitudinal 
wave from the satellite waves, which are shap e d formed in the time of the reflection of last from 
free earth surface. 

The recording of the displacement vectors of the direct longitudinal (or shear) wave is are 
carried out simultaneously in the investigated and check boreholes, moreover in the investigated 
borehole - with the help of the three-component multimodular downhole multi-level borehole sonde 
sequentially moving along the depth of the borehole depth , and in the check borehole with th e h e lp 
by means of the stationary located three-component one modular level sonde. The detailed three- 
component profiling in the investigated borehole is executed with the step of the discrete 
observations every (2.5 - 5) m. 

Since the s e ismic wav e s propagat e outsid e th e boundari e s of th e inv e stigat e d bor e hol e and 
insonity the vicinity of the borehol e , th e r e liabl e diagnostics of such capacitativ e , filtrational and 
fluid saturated prop e rti e s of the r e s e rvoirs in th e conditions of th e ir natural b e dding is carried out 
monotype time signals extraction for the compression and shear constituents out of the observed 
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experimentally seismic wavelets on the axial and tangential components of the displacements vector 
for the direct longitudinal wave is realized - for the compression constituent by means of deduction 
the seismic wavelets, observed on the axial component with weight coefficient out of the seismic 
wavelets, observed on the axial component whereas for receipt of the shear constituent - seismic 
wavelets, observed on the tangential component are normed with the same weight coefficient for 
each observation point in the investigated and check boreholes . 

The interval values calculation of the elasticity modules (of longitudinal and shear) is made by 
means of the coprocessing of longitudinal and shear waves velocities and the density logging data. 
The decoding is realized by means of the parametrical analysis of the monotype time signals for the 
compression and shear constituents of the direct longitudinal wave with receipt of the quantitative 
estimates of their significant dynamical parameters for the borehole and check data, moreover in 
quality of the significant dynamical parameters of the monotype time signals for the compression 
and shear constituents of the direct longitudinal wave makes use of the initial amplitudes, steepness 
exponents, attenuation coefficients, instantaneous frequencies and initial phases. 

With the purpose of the receipt of dynamical parameters of the layers impulse responses for 
given point of the observation executes numerical correction of the dynamical parameters estimates 
of the monotype time signals for borehole observations by means of the utilization of the analogous 
parameters of signals for the check observations, and further by means of the layer-stripping 
method is realized the correction of revised estimates of the dynamical parameters for the given 
observation point by means of use the analogous parameters of the monotype time signals for the 
previous observation point of lengthways depth of the borehole. For qualitative indication of the 
presence and fluid type - oil, gas or water, filling the pore and fissured volume in the rocks, 
quantitative values of the bedded absorption coefficients for the compression constituent of the 
direct longitudinal wave are compared with the calibrated values of the above-mentioned 
coefficients for the oil-gas-or water-saturated rocks. 
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The calculation of bedded coefficients of the effective attenuation is realized by means of 
statistical processing of the interval values of the initial amplitudes previously corrected owing to 
the conditions change in the time of multiple excitation of the oscillations, and owing to the 
geometrical divergence of the wavefront, further is defined the bedded attenuation coefficients by 
means of the statistical processing of the interval values of the longitudinal modules of the 
elasticity, and also is defined bedded absorption coefficients, for that the absorption data average 
out statistically on depths intervals of the layers, accordingly for the compression and shear 
constituents of the direct longitudinal wave. For the calculation of the oil-field parameters collection 
of the reservoirs, such as bedded values of the coefficients of porosity, granular and fissured 
permeabilities, residual water saturation and fluid saturation, are used bedded values of the 
coefficients of volumetric and shear dynamical viscosity jointly with values of the temperatures, 
diffusion constants, meanderingness of pore channels coefficients and others. 

Under the given invention all the significant parameters of the reservoir and fluid-saturated 
properties of the rocks in the radius of £50 - 100} m for the vicinity of the borehole are received, 
whereas the above-mentioned parameters are not defined with the required completeness by the 
standard complex of the well logging methods, which give only the restricted information about the 
vicinity of the borehole in the radius which is not exceeding (0.5 - 2.5) m. 

With the help of the given invention the reliable information about the quantitative 
parameters of the oil-and-gas-bearing layers can be received, irrespective of the difficulties of the 
drilling and bedding conditions, the lithological content and the thermodynamical state of the oil- 
and-gas deposits, that finally reduces the cost price of the drilling expenses and owing to the cutting 
of the common number of the boreholes maintains the ecological situation in the region of 
exploration. 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
For the more complete understanding of the main point of the given invention and its 

advantages the references to the description are made, which is considered here together with the 

accompanying figures where the following is represented: 

Figure 1 - the seismograms , corrected across beyond the first arrivals, of the compression 

(S) and shear /T) constituents of the direct longitudinal wave for th e compr e ssion (S) and 

sh e ar (T) compon e nts according to th e data of th e downhol e obs e rvation b orehole (A) and check 

ob s ervation (B) of 3C VSP for observations ( the borehole of the Lopushna-1 3). 

Figure 2 - the calculated dynamical parameters of the impulse responses of the rocks for the 
borehole of the Lopushna-1 3 : the attenuation coefficients of t e mporal att e nuation compression f fi Q ) 

and coefficients of shear attenuation fa ): the frequencies of natural fr e qu e ncy oscillations (F 0 ), 

which are compared te with the lithology. 

Figure 3 - the calculated rheological parameters of the rocks and parameters of the initial 
amplitudes effective attenuation of seismic wavelets for the borehole of the Lopushna-1 3 : th e e lastic 

moduli (Lam e 's contans) — of the compression (A)and shear (^i); t he coefficients of volumetric ( X ) and 
shear (tj) dynamical viscosity jointly with the initial amplitudes diagrams of \nA' c f, fHj] b\ means of 

which defines bedded coefficients of the effective attenuation fa'" att (lnA' c '" which are compared 

te -with the lithology. 

Figure 4 - the calculat e d att e nuation and absorption parameters and the elasticity modules of 
the rocks for the borehole of the Lopushna-1 3: fa - the absorption coefficients (for the 

characteristic frequency c5 - 27c20(Hz ) -a'dflA' p) and diagrams of the elastic modules 
logarithms Ln( X + 2jj, )[H] a Ln/u(H ) by means of which the bedded attenuation coefficients (a'- " ) 
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are defined, characterizing gradient properties, which are compared te-with the lithology. The 
param e t e r D a sampling varianc e ; 

Figure 5 - the calculat e d reservoir and fluid saturated saturation parameters of the rocks for the 

borehole of the Lopushna-13 : the coefficients of porosity ( ® g ) 9 granular permeability (K g ), residual 
water saturation (£ ! rvv )and fluid saturation (k' f s ) with the d e finition of th e fluid z type 

identification , filling the pore spac e volume of the reservoir, which are compared to the lithology. 

Figure 6 - the calculat e d reservoir and fluid saturat e d saturation parameters of the rocks for the 
borehole of the Lopushna-13 : the average length of fractures (I), the linear density of fractures ratio 

(L df ), the coefficients of fissured permeability (K f ) and fluid saturation (k" fs ) with the 

d e finition of the fluid z type identification, filling the fissured spac e volume of the reservoir, which 
are compared to the lithology. 

DETAILED DESCRIPTION OF THE INVENTION 

The diagnostics of the dynamical, rheological, absorbing, reservoir and fluid-saturated 
properties for the rocks in the vicinity of the borehole is based on the combined deformation model of 
the vertically - inhomogeneous viscoelastic isotropic solid medium for the isothermal case, proposed 
by the author [Briginevich V.A. 1991. The combined model of deformation and the conditions of 
longitudinal wave propagation in a viscoelastic solid. Geophys. Journal. 10(3), 388-405]. The 
generalization of the above-mentioned deformation model of the viscoelastic solid medium for case 
of the verticallv-inhomogeneous changing of its parameters over the depth are discussed below. 

According to the combined model at th e d e forming of deformation model of the vertically- 
inhomogeneous viscoelastic solid medium such obs e rvationally m e asur e d macroscopical macroscopic 
parameters as the displacements, stresses, e lastistic moduli modules of elasticity , the coefficients of 
dynamical viscosity and others dir e ctly dep e nd on th e macroscopical can be measured. Moreover above- 
mentioned macroscopic parameters directly depend on such microscopical parameters of the real solid 
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medium as microstructure of the matrix material (grains) , the sizes of the micropores, the density of 
microfissures and their fluid saturation, mor e ov e r thoir activity is shown through th e when their presence 
is appeared by means of total response of the last to th e e nclos e d e xt e rior impulse external action in 
the shap e form o f the average statistical values of the-porosity, Assuring, permeability and fluid 
saturations saturation . 

The combined deformation model of the vertically-inhomogeneous viscoelastic solid 
medium is determined by the equations, which for the rectangular coordinate system (x,y,z) 

characterize the following stressed and deformed state: a. >a x =<r y ; a xz =0^ =a zx = 0; £ : *Q; 

e x =e y = 0 9 where a x y z - normal stresses; , , - tangential stresses; e x y z - linear 

deformations, which arise in the solid medium during the propagation of plane longitudinal wave, 
moreover the equations are represented for the components of the stress vector in the operator form: 



°z = 



X + X — +2/i 



dt 



dt, 



d^ 
— + — 



U dt 



J 



(1) 



A + A. — 
dt) 



(2) 



where X = X(z), fi- /u(z) 



the compression and shear moduli (Lame's 



coefficients); X* = X*(z), tj = t](z) - the coefficients of the volumetric and shear dynamical viscosity, 
respectively. According to the form of the e quations of th e state equations (1,2), which satisfy the 

movement equation: ^^ = p^-^-, where p = p(z) - the density; 1 = (x,y,z) we receive the 

dx t dr 

homogeneous wave equations for the displacements components of the plane longitudinal wave 
propagating in the direction of the axisoz : 



H dt 



f d 1 „ d^ 

— + 2a — 
{dt 2 dtj, 



_ + (2« + /0- + « 



>U. = V: 



1 2 + ajn(\l + 2/i) d 



dz 



dz 



dz 



(3) 



P — r+ A 
dt 1 



dz 2 dt 



+ x 



dz 



U = 0 , 

x,y ' 



(4) 
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where u =U w (z,t) - the components of the displacements vector for rectangular 
axe&(x,y,z); V = V p (z) = [(Z + 2ju)/ p] I/2 - the local velocity of the longitudinal wave; 



a> o =(Xp/ £ri) in - the angular frequency of natural oscillations; a ~ V / 2r l - the shear coefficient of 

attenuation; p = (X + 2^i)/ X - the volumetric coefficient of attenuation. 

Within the the correct statement of th e correctly formulated nonstacionary mixed problem 
on radiation in case of the inhomogeneous initial and homogeneous boundary conditions of the 
impedance type which are executed on the flat plane wavefront of the direct longitudinal wave, 
when t > 0, z > 0, the solution of the equations (3,4) or the impulse responses of the deformation 
for the above - mentioned medium, using the m e thod of Fourier method, have been received for the 



displacements component of displac e m e nts in the shape: 



U z (z,t)= A ° 



(I + 2 m /' 



— ■ - B 0 e' a sin(cot + 
Pa 



U, Jz.t) = -A 0 • B 0 e' m ■ sin( at + <p 0 )- exp(-k k z). 



(6) 



exp( ikz ), 



where A 0 = const ; a> = (a> 0 2 -a 2 ) 1 ' 2 - the angular frequency of attenuaton oscillations; 
P 0 =(o 0 2 /p- the coefficient of volumetric attenuation; B 0 =(a 2 +co 2 ) 1/2 /ao) - the amplitudes; 



<p 0 = arctg(o) / a) - the initial phase; k = k(z) = j 



dz 



dz; k p =k p (z) = G) 0 /V p 



z 

the local wavenumber; k k = \k x dz\ k x = k x (z) = co 0 /V x - the absorption coefficient; 

0 

V x =V x (z) = (X/ p ) 1/2 ; i - the imaginary unit. 

The above-mentioned solution of a direct problem allows to predict the superposition of two 
monotype t e mporal time signals in the direction of the propagation of the direct flat plane 
longitudinal wave for each point of the vertically-inhomogeneous viscoelastic solid medium. The 
presence of the dual impulse response during the longitudinal wave propagation is defined by the 
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conditions of the stress-deformed state of the liquid-porous (for compression constituent) and 
fracture-brittle (for shear constituent) elements of the dynamically deformable solid medium. 

Therefore, the monotype time signals which are excited in the continuous solid medium by 
the flat plane longitudinal wave for each point of th e vertically inhomog e neous viscoelastic solid 
m e dium. Th e presenc e s of th e dual impuls e response during th e longitudinal wav e propagation is 
d e fin e d by th e conditions of th e str e ss d e formed stat e of th e liquid porous (for compression 
compon e nt) and fractur e brittle (for shear component) e l e m e nts of th e dynamically deformabl e 
solid m e dium. Th e r e for e , the monotyp e temporal signals which arc excited in th e continuous solid 
medium by th e flat longitudinal wav e, will consist of two damped processes differing by the 
polarities, initial amplitudes, t e mporal attenuation coefficients of att e nuation , instantaneous 
frequencies and initial phases, and only one of these processes can be clearly observed on the shear 
components of the displacement vector. Really, from the experimental data of the VSP polarization 
method of VSP is known, that the direct longitudinal wave propagating in the stratified solid 
medium is polarized linearly for the first and is nonlinearly for the subsequent phases of the seismic 
impuls e s wavelets (Galperin E.I. 1984. Polarization method of seismic prospecting. D. Reidel 
Publications Co). 

It is provided that on the flat plane wavefront of the longitudinal wave propagating in the 
direction of the axis oz , the following boundary coditions of the stationary type for the component 



U z (z>t) ar e carried out: 



d 2 ^dln(A, + 2fii) d 



V 



dz 



dz 



+ k 2 
dz kp 



J 



UJz) 



< =0, the equation (3) is converted to 



the form: 



d 2 ^dln(k 



dz 2 



dz dz j^dt 2 



a p 



d 2 dln(X + 2fx)d^ 



V 



dz 2 



dz dz 



U z (z,t)=0. 
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From this equation according to the Bodgia theorem we received: 



H dt p 



dln(X + 2n) d 



V 



dz 



dz 



dz 



J 



U' s (z,t)=0;(7) 



dt 2 dt p 



f d 2 i dln(A + 2y) d 



V 



dz- 



dz 



dz 



J 



u;{z,t)=o,(8) 



where U 2 {z,t)=U' z (z,t)+u:(z,t). 

If the boundary conditions are executed, each of the equations (7), (8) will fulfil to the 
sample solutions for the harmonic components: U z ' r "(z,t) = e ~ lcot exj^k(aj)z\ (9), where the 
characteristic angular frequency c5 will be the independent variable; k{c5) - the complex 
wavenumber; k(oJ) = co IV(c5) - the real wavenumber; a (co) - the attenuation coefficient; K(<y)-the 
propagation velocities of the harmonic components. Taking into account the second approximation 
of the VKBD method (*(*>) = k p ) 9 the solutions (9) will be true if: V{W) = V p m /<a 09 a u "{a)=a' "{c5) 

+ a'''\ where V(cd),a' "(aj) - the propagation velocity and the total attenuation coefficients, 



respectively; a'(a>) = Pco /2V o) 0 and a"(co ) = a(Q / V co 0 



the absorption coefficients, 



characterizing the rheological (visco e lastic) properties of the solid medium; a ft "=dln( X + 2jj)/ dz 

- the attenuation coefficient, characterizing gradient properties of the elastic parameters in the solid 
medium for U' 2 (z,t) and U" z (z,/), respectively. 

The description of the propagation of the direct flat plane longitudinal wave in each 
observation point of the stratified vertically-inhomogeneous viscoelastic solid medium is possible 
while using the convolutional model and approximations approximation of the seismi c wavelets by 
means of the analytical function of Berlage representing the description of the wide class of the 
monotype t e mporal time signals. 

Hence the impulse responses of the above-mentioned deformation medium for the 



components of the displacement vector may be offered in the form: 
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u x y (z,t) = -A } " t exp(-at ) sin(aj" t - (p 0 ")exp{^ k x z} 

< 

u z (z,t) = \ij 't y exp(- fi 0 t)sin(o) f t -p 0 ')- A } "t p " exp(- at)sin(co" t - (p 0 ")]exp(ikz), 
where u xyz (z,t) - the displacements for the (x,y,z) - components; A/, A/' - the amplitudes 

dependent on the coefficients of the wave transmission in the layers; p' " - the e xponents of 
steepness exponents ; p o ,a - the t e mporal co e ff e ici e ntG of attenuation coefficients ; a' "- the 
angular frequencies; ^/"- the initial phases for the compression (') and shear fl constituents, 
respectively. 

The physical observations of the three-component VSP data for the recording of the 
displacement vector are realized in the vertical cylindrical boreholes which are located in the stratified 
vertically-inhomogeneous viscoelastic solid medium. According to the classical mechanics of the 
continuous medium the borehole is the concentrator of stresses and deformations in the above- 
mentioned medium. It is known, that in the presence of orthogonality of the axis oz the boreholes 
concerning the wavefront of the flat plane longitudinal wave, the wavelength considerably exceeds 
diameter of the borehole 2a, where a - the radius of the borehole. Therefore, for the wavefield 
components and observations on the wall of the cylindrical borehole the author has received the 
special solution of the problem, where in the calculation is based on the following. 

If the vicinity of the borehole is considered in the cylindrical coordinate system with the 
ordinate oz , directional along an axis of the borehole, then the boundary conditions on the borehole 
wall for the stresses are usually set in the form: a r \ r=a = 0 ;a e \ r=a = 2<r 0 ;cr z \ r=a = a 2 , where <j r ,<j 0 ,cr z 

- the normal stresses for the radial, tangential and axial components. According to the received 
solution the components of the displacement vector which are formed by the direct flat plane 
longitudinal wave on the borehole wall, propagating in the vicinity of the borehole in the direction 
of the axis oz, are: 
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u r (z,t ) = 0; u 0 (z,t ) = -j2Aj"t r exp(-at )sin{a>"t- (p 0 ")exp{- k x z); 

u z (z,t) = exp(- fi 0 t)sin((o't'-<p 0 f )-A J "t r exp(-at)sin(3"t -y 0 "^exp{ikz\ 

where the ind e x (~) u re /z.t ) - the displacements components in the cylindrical 
coordinate system (r,0,z) ; the index (~) designates the estimated parameters during processing . 

Then at excitation of the direct longitudinal wave in the vicinity of the borehole and during 
the displacements observations on the borehole wall the superposition of the monotype t e mporal 
time signals on the axial component is defined: u z (z,t)\ z=H =u z (t) = u z '(t) + u z "(t) , whereas on the 

tangential (or transverse) component: u e (z, t )\ Z=H =u e ( t ) only one monotype temporal time signal is 

observed: u T (t) = u z (t)-u e (t)/^j~2 . Hence, it follows: u s (t) = u z (t)-u 0 (t)/ 42 , where u z (t),u e (t) _z 

the displacements observed on the axial and tangential components; u s (t),u T (t) - the 

displacements characterizing the monotype t e mporal time signals for the compression (S) and 

shear (T) compon e nts, accordingly constituents of the direct longitudinal wave, respectively . 

The process of the utilization of th e invention utilization includes: the observations of the 
wavefield in investigated and check boreholes with necessary fullness ( r e gistration of the 
displac e m e nt v e ctors by means of 3C VSP three-component observations for recording the 
displacement vector) and with necessary detail (the compact dense step for Ae-observation points at 
th e profiling of th e inv e stigat e d during the borehole profiling) ; the digital recording of seismic 
waves in the at wide-band dynamical dynamic and frequency ranges; the data processing of 3C VSP 
three-component observations with the purpose of their transformation from the obs e rvation point 
ef view te of the data paused, specified by the-observation system (the near surfac e subsurface 
source of the for wave excitation, while the-geophones are in the borehole) A to the aspect of the-data 
which represent values of the dynamical parameters of the impulse responses of the rocks rock 
layers , and th e n will b e carried out by the g e ological subsequent geologic interpretation by means 
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of the p e trophyscial rheological data converting to the values of the p e trol e um oil-field parameters 
for th e lay e rs of the rocks. 

The procedures of the downhol e borehole observations envisage the recording of the three- 
component data in two boreholes simultaneously, moreover one of them - the investigated borehole 
where measurements are realized by means of profiling of the l e ngthways d e pth of th e b orehole 
depths intervals with the three-component multimodular downhol e multi-level borehole sonde, and 
another - the check borehole in which the three-component one modular level sonde is located 
stationary. 

For the oscillation excitation one near-surface source (the shotpoint) is used, moreover the 
distance of the last from the investigated borehole does not exceed (50-100) m. The oscillation 
excitation is realized by the known method, for example, by means of the location in the near-surface 
borehole of the explosive charge (or the air gun) at the optimum depth which guarantees the 
receiption reception of the maximal energy and the simple shape of the seismic wavelets of the direct 
longitudinal wave. 

The optimum depth of the location of the explosive charge (or the air gun) is determined by 
means of the execution of the obligatory condition of the separation of the direct longitudinal wave 
from the satellite waves, which are shaped in the time of reflection from the free surface. 
Simultaneously the force of the impulsive action at the excitation of oscillations, determined by the 
weight of the charge (or the air gun power) and the seismogeological conditions of its location, 
should guarantee the reception of such energy of the direct longitudinal wave in the investigated 
borehole in order to enable the seismogram of the tangential component for the downhol e borehole 
three-component sonde to distinguish the useful signal in the background of the noise (the signal-to- 
noise ratio must exceed 1). 

For the registration of the conditions change of the nonsingle excitation of oscillations 
during the production of works on the investigated borehole the check borehole is used, located 
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close to the area of the excitation and on the trajectory of the direct waves between the investigated 
borehole and the shotpoint at the distance of (20-50) m from th e latt e r last . The three-component 
one modular level sonde is dipped in the check borehole to the depth below the depths of the 
location of the explosive charge (or the air gun). 

The borehole investigations are conducted by means of the standard equipment with the 
digital recording of the seismic signals, which implements on the open channel with the sampling 
interval of 0.5 ms and in the time of absence of the amplitude adjustments. The investigations are 
conducted in the borehole at the perspective in the oil-and-gas-bearing relation interval of depths 
(length 300 - 500 m) by means of the detailed three-component profiling the latt e r last with the step 
of the discretic observations every (2.5 - 5} m. 

The method of the data processing of the downhole borehole observations is based on the 
solution of the inverse problem for the case of normal plane waves incidence. The method envisages the 
data processing of the three-component observations which preliminarily are oriented in the space, and 
the extraction of the monotype t e mporal time signals for the compression and shear constituents from 
the seismic wavelets which are observed experimentally on the axial and tangential components of the 
displacement vector of the direct longitudinal wave for the downhol e borehole and check devices. 
The procedures of the digital extraction of the compon e nts constituents and respective monotype 
t e mporal time signals were considered earlier. 

On the first stage of processing by means of the parametric analysis of the seismic traces 
(Taner M.T., Koehler F. and Sheriff R.E. 1979. Complex seismic trace analysis. Geophysics \ 44, 
1041 - 1063.) the monotype temporal time signals are decoded according to the compression and 
shear compon e nts constituents of the longitudinal wave with the further receiving of the 
quantitative estimates of their significant dynamical parameters such as the initial amplitud e 
amplitudes , th e e xpon e nt of steepness exponents , th e t e mporal co ee fici e nts of attenuation 
coefficients , the-instantaneous fr e qu e ncy frequencies and initial phas e phase s. 
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The parametric analysis of the monotype t e mporal time signals tog e th e r jointly with the 
procedure of the determination of the numerical estimates of their dynamical parameters is 
performed with the use of the approximation of the seismic wavelets by the analytic B e rlag e 
function of Berlage which allows to represent adequately the typical monotype t e mporal time 
signals by means of the totality of only five significant dynamical parameters: 

U SJ (zj) = A 0 '"t exp[- (p 0 , S)t]sm(2aF" "t + $ 0 "") 9 
where the estimated parameters: A 0 ''" - the initial amplitudes; p' " - the e xpon e nt of 
steepness exponent ; f3 0 ,a - the t e mporal co ee fici e nts of attenuation coefficients ; F r " - the 

instantaneous frequencies and <p 0 ' " - the initial phases of the monotype t e mporal time signals. 

In the invention the opportunity of receiving of the most precise values of the parameter is 
taken into account, since each of the dynamical parameters at measuring its value will have biased 
estimates. The latt e r last is achieved by means of the computer modelling of the monotype t e mporal 
time signals, according to which the received quantitative estimated parameters by means of the 
model of the synthetic signals are compared visually to the observed monotype t e mporal time signal 
and the next are the corrected with the application of the optimal deciding rules. 

It is important, that the prognosis of the shape of the monotype temporal time signals 
ensures the receiving of the estimates of their dynamical parameters for each observation point of 
the lengthways depth of the borehole. The estimated parameters can be numerically corrected for 
converting them from the observable values into the typ e of bedded values type . For this purpose, 
using the obs e rvation downhol e borehole and check data of the observation for the respective 
component of the direct longitudinal th e obs e rvation downhol e and check data for th e r e spectiv e 
compon e nt constituent of the direct longitudinal wave, they separately calculate the estimated 
dynamical parameters of the monotype t e mporal time signals, and after the improvement of the 
received values of the estimates by means of the computer modelling execute the correction of their 
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values according to the algorithms, offered by the author, by the recalculation of the parameters, 
characterizing the impulse responses of the rock layers. 

G aining th e dynamical parameters of th e monotyp e t e mporal signals according to th e 
postulat e of th e lin earity-In view of the postulated linearity of conditions for displacements (in the 
time of the seismic waves propagation), prop e rti e s of r e civing receiving and recording devices^ the 
problem of the linear filtration dynamical parameters acquisition of the monotype time signals can 
be considere d, mor e ov e r within the limits theory of the linear filtration, when the seismic channel 
has the impulse and frequency characterizations characteristics . For exception of the distorting 
effect owing to the activity of the recording channel and for the exception of the conditions change 
of the wave excitation at the beginning the numerical correction of the estimated dynamical 

parameters ( P 0 ,a;Sf" ) for the monotype t e mporal time signals of the downhol e b orehole 

observations is implemented by means of the use of the analogical parameters out of the check 
observation data according to the dependencies: 

( Packer = 2 ( Po> S )bor ~ ( ' Po>% kheck' <°lor = 2 &lor ' Kheck ' 

where (fi 0 ,a ;5" )bor. and (]} 0t a ;S" ) C heck- the estimated values of the dynamical parameters 

for the downhole borehole and check data and the respective observation point. Then the correction 
of the revised estimates of the dynamical parameters is implemented by means of the layer-stripping 

method ( P 0 ,a;G> n J^, that allows to delete the distortions of the parameters for the downhol e 
borehole data, caused by the filtrational effect of the covering thickness of the rocks, and also to 
convert the latt e r last p arameters into the respective parameters of the impuls e r e spons e s of layers 
impulse responses . 

The recalculation of the revised estimates in the interval values of the dynamical parameters 
of the impulse response for the separate layer is carried out according to the following formulas: 
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where (P 0 ,a; p' " ,g>" )j - the values of the revised estimates of the dynamical parameters for 
the j -observation point; ( f} o ,a;p''" t cf ) M - the revised estimates values of the dynamical 
parameters for the (j-1) - observation point lengthways depth of the borehole; (P 0 ,a t (Q n ) layer - the 

interval values of the dynamical parameters of the impulse response of the layer. 

On the second stag e of th e processing th e int e rval valu e s of th e e lastic moduli ar e d e fin e d 
which ar e calculat e d by m e ans of th e use of th e kin e matic param e t e rs of th e longitudinal and shear (or 
conv e rted) waves. The joint use of the kinematic parameters and th e d e nsity log of longitudinal and 
shear waves, and also the density logging data allows to d e fin e th e num e rical on the second stage of 

the processing allows to define interval values of such elastic moduli modules as longitudinal - 

■j 

(X + 2fj,) = pV p ~ an( j s hear -p = pV s 2 , where A,//- Lame's coefficients; V pj - the velocities of 

longitudinal and shear waves and P - the density. Using the values of the dynamical parameters of the 
impulse responses and the elastic moduli modules for the rock layers, the values of such rheological 
parameters as the coefficients of the dynamical viscosity: X = (k + 2fi)/ p (for the volumetric viscosity) 

and ? = M ' 2a (for the shear viscosity) are calculated. 

On th e third stag e for th e pr e viously s e parat e d lay e rs (according to th e lithologic data) th e 
sp e cial group proc e ssing of th e e stimat e d valu e s of initial amplitud e s (.sub.OV) and th e longitudinal 
e lastic modulus (.lambda.+2.mu.) ar e impl e m e nted. Th e giv e n invention e nvisag e s the possibiliti e s 
for incr e asing th e accuracy of th e obtaining of th e amplitude e stimat e s of th e s e ismic impuls e s of 
th e direct longitudinal wave: by means of th e us e of th e monotype t e mporal signals; obtaining their 
estimates of th e initial amplitud e s (.sub.O'"); th e corr e ctions of e stimat e d valu e s of th e abov e 
mention e d amplitud e s for th e d e cr e as e of th e ir straggling owing to th e chang e of th e conditions in 
the time of th e multipl e wav e e xcitation and th e ir subs e qu e nt normalization for th e comp e nsation of 
the geom e trical div e rgenc e of the w r av e front (owing to th e differ e nc e of th e latt e r from ideally flat) 
for th e compr e ssion and sh e ar components of th e direct longitudinal wav e , accordingly. 
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[0070] Aft e r th e calculation of tho calibration coefficients out of th e check obs e rvation data and the 
correction of the estimated values of th e initial amplitudes for th e downhol e obs e rvation data th e 
normalization of th e above mention e d amplitudes is performed with th e aim of the comp e nsation of 
th e g e om e trical div e rgenc e of th e wav e front of tho dir e ct longitudinal wav e . Furth e r, according to 
th e revis e d and normaliz e d valu e s of th e initial amplitud e s for th e e ach obs e rvation point 
l e ngthways d e pth of th e borehol e and for th e longitudinal wav e the resp e ctiv e component the 
co e ffici e nt of th e eff e ctiv e att e nuations for th e j layer io calculated: V'.sub. e ffant.~(ln A.sub.0j2'," 
In A.sub.OjlV^DELTA.H.sub.j, wh e r e A.sub.0jl m , A.sub.OjT" th o valu e s of amplitudes on th e 
upper (j.sub.l) and low e r (j.sub.2) boundari e s of th e j lay e r with th e thickn e ss .DELTA.H.sub.j, 
which ar e d e fin e d statistically by m e ans of th e l e ast squar e s m e thod using th e d e clination of th e 
profil e of th e chang e of th e natural logarithm of amplitud e s (inA.sub.Ojl 1 ,") for th e e stimated 
int e rval e£ d e pths. 

[007 1 ] Th e co e ffici e nts of att e nuation for th e dir e ct longitudinal wav e according to th e combin e d 
d e formation mod e l of th e v e rtically inhomogen e ous visco e lastic solid m e dium ar e as follows: (tild e 
ov e r (.alpha.)) '/'((ov e rscor e (.om e ga.)) )- (ov e rscor e (.alpha.)) , ,"+.alpha. f / , ( (ov e rscor e (.om e ga.)) ), 
where (ov e rscor e (.alpha.))'," th e att e nuation co e ffici e nts characterizing th e gradi e nt properti e s of 
th e e lastic param e ters; .alpha.',"( [ov e rscor e (.om e ga.))) th e absorption co e ffici e nts charact e rizing 
th e rheological prop e rti e s of th e lay e rs of th e rocks for th e compression and sh e ar compon e nts, 
respectiv e ly. 

[0072] The real solid m e dium usually is thin lay e ring which at th e d e finite assumption for th e 
s e parat e thick lay e r can b e regard e d as the gradient (or v e rtically inhomog e n e ous) m e dium. Th e n 
the co e ffici e nt of att e nuation .arpha.',"( (ov e rscor e (.om e ga.))), charact e rizing th e e ffects of e n e rgy 
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disp e rsion in thin layering b e ds, can b e calculated by moans of th e l e ast squar e s m e thod according 

tei (ovorscor e (.alpha,)) , ,"^ln(.lambd a.+2.mu.).subj2 

ln(.lambda.+2.mu.).subjl]/.DELTA.H.sub.j, — wher e — (.lambda.+2,mu.).sub.j the valu e s — of th e 
longitudinal elastic modulu s on th e upper (j.sub.l) and lower (j.sub.2) boundari es of th e j layer with 
tho thickn e ss .DELTA.H.sub.j. Th e distinctiv e f e atur e of th e m e thod for th e calculation of th e 
absorption co e fficients .alpha.V'(cs) in th e giv e n inv e ntion th e r e is th e opportunity of th e us e of 
only valu e s of th e t e mporal dynamical param e t e rs (.beta., .alpha., .om e ga..sub.o) and th e velociti e s 
of th e longitudinal wav e Vp in th e lay e r to th e obtaining of th e above mention e d absorption 
co e ffici e nts. 

[0073] Th e absorption co e ffici e nts, lin e arly d e p e nd e nt from th e charact e ristic angular fr e quency 
(ov e rscor e (.om e ga.)}, ar e calculat e d for e ach obs e rvation point l e ngthways depth of tho bor e hol e 

according te formulas: .alpha.'( (overscor e (.om e ga.)) )~.b e ta. (ov e rscore 

(.om e ga.)) /2V.sub.p.om e ga..sub.o; '"( (ov e rscor e (.om e ga.)) ))- alpha, (ov e r score 

(.om e ga.)) /V.sub.p.om e ga..sub.o, and th e quality factors (Q factors) ar e calculat e d according to: 

Q'( (ov e rscor e (.om e ga.)} )~.om e ga..sub.o a2/.b e ta. (ov e rscor e (.b e ta.)}; Q"( (ov e rscor e 

(.omega.)) ~.om e ga..sub.o/2.a — lpha. (ov e rscor e — (.om e ga.)) — for — the — compr e ssion — and — sh e ar 
compon e nts, accordingly, wh e r e .omega.. sub.o~2 .pi.F.sub.o th e angular fr e qu e ncy of natural 
oscillations; .b e ta .-.om e ga.. sub.o.sup.2/D. th e t e mporal co e fficient of volum e tric attenuation. Th e 
rec e iv e d int e rval valu e s of th e absorption co e ffici e nts .alpha. ',"( (ov e rscor e (.omega.))) lat e r ar e 
av e rag e d statistically by m e ans of th e least squar e s m e thod on th e int e rval of th e d e pths of th e 
layers. 

Simultaneously the absorption coefficients, linearly-dependent from the characteristic angular 
frequency ( a) ), are calculated for each observation point lengthways depth of the borehole 
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according to formulas: a'Cco ) = p^/2V p a> Q ; a"{aj) = ac5/V p G) 0 , and the quality factors ( O -factors) 

are calculated according to: OYco ) = (o 0 2 / Po5;Q n (o5) = a> 0 2 /2aco , respectively for the compression 
and shear constituents, where c5 < (o c - the characteristic angular frequency; co c = 2nF c - the angular 
frequency of natural oscillations; 0 = co 2 / B c - the volumetric attenuation coefficient. In addition 
distinctiveness of the calculation method of the absorption coefficients \ a'- "(co) ] in the given 
invention is utilization for their receipt only dynamical parameters ( 0.a.a) c ) and the velocities 
values (V p ) of longitudinal wave. 

The author has established experimentally that the bedded values of the absorption 
coefficient Lla f (dJ)\ for the compression component constituent of the longitudinal wave can serve 

as the indicator of the hydrocarbon presence in the r e servoir, sinc e th e s e ismic non e lasticity of th e 
last is condition e d by such its rocks what is caused their properties such as the porosity, fracturing, 
permeability and fluid saturation (Brygynevych V. A. 1997. REOKONA seismic technology studies 
reservoir and fluid - saturated properties of the rocks in the vicinity of borehole. 59 EAGE 
Conference, Geneva, Switzerland, Extended Abstracts, P075). Therefore, the calibrat e d bedded 
values of the absorption coefficient {a'((o)\ for the compression compon e nt constituent of the 

longitudinal wave are used to receive the quality estimates of the fluid z type identification, filling the 
pore and fissure volume in the rocks. For this purpose quantitative values of the obtained absorption 
coefficients for the compression constituent reduced to bedded conditions, are compared with the 
calibrated values of the above-mentioned coefficients for the oil-gas-or-water-saturated rocks. 

In this connection the calibration conditions are offered: a'Jco )<a' 0 (<5 )<a' g (a) ) , where the 

quantitative values of the absorption coefficient are compared: for water - a'Jaj ) , oil - a r 0 (o5 ) and gas - 
a'(W), respectively. The calibrated value of the absorption coefficient for the compression constituents 
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depends from the characteristic frequency (co). Thus, for the frequency m =2n20(Hz) the numerical 

3 1 * 

values of the absorption coefficient changes within limits: for water - a' w (cJ)= (1-4)*10* m" , for oil - 
a' 0 (o5) = (5-7)*10" 3 m _1 and for gas - a' g (W) = (8 - lO^lO^m' 1 . 

In th e ir tu m^-However, the calibrated values of the absorption coefficient [a n (cd)\ for the 

shear compon e nt constituent of the longitudinal wave are changed in the reverse order: 
a"Ja> ) > a n 0 (c5 )>a" g (oj) , where the absorption coefficients are: for water - a n Jo5 ^ , oil - a n 0 (aj ) and 

gas - a"Jco), respectively. It has been determined by the changes of the fluid type, filling the por e s 

pore and fractures fissure volume in the rocks, mor e ov e r the absorption coefficients [c5'' "(co )\ for 

the compression and shear compon e nts chang e constituents are changed in the opposite directions. 

The given invention envisages an opportunity of the precision increase of receipt of the 
interval initial amplitudes estimates for seismic wavelets of the direct longitudinal wave by means 

of: the monotype time signals utilizations; presence of the initial amplitudes estimates ( A^") ; the 

corrections of the above-mentioned amplitudes estimates for decrease of their scatter owing to the 
change of the conditions in the time of multiple excitation of the oscillations and their subsequent 
normalization for the compensation of geometrical divergence of the wavefront (owing to 
difference of last from plane). 

On the third stage of processing after correction factors calculation out of the check 
observations data and the subsequent correction of the estimates values of the initial amplitudes for 
the borehole observations data is carried out normalization of revised amplitudes for the 
compensation of geometrical divergence of the wavefront respectively for the compression and 
shear constituents of the direct longitudinal wave and the corresponding observation point. 
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Then for the previously separated layers (according to the lithological data) the statistical 
processing of the interval values of the initial amplitudes (A' 0 ") 3 the longitudinal modules of the 

elasticity (X + 2u) and the absorption coefficients \a' "(co )] is implemented. 

Further, according to the revised and normalized values of the initial amplitudes for the each 
observation point lengthways depth of the borehole and for respective constituent of longitudinal 
wave the effective attenuations coefficient for the /'- layer is calculated: 

a'-" effatt =(lnA 0 '" J2 -lnA c '" jj )/AH j % where A c ''" jn A C ''" J2 - the values of amplitudes on the 

u pper f / ; j and lower ( j 2 ) boundaries of the / - layer with the thickness (AH, ) , which are defined 

statistically by means of the least squares method using the declination of the changes diagrams of 
the natural logarithm of the initial amplitudes ftnA c ' " j ) for the estimated depths interval. 

The total attenuation coefficients for the direct longitudinal wave according to the combined 
deformation model of the vertically-inhomogeneous viscoelastic solid medium are follows: 
a' "(c5 ) = a' " + a' "((D) 3 where a' " - the attenuation coefficients, characterizing gradient properties of 

the elastic parameters; a f ' "(75 J - the absorption coefficients, characterizing Theological properties of 

the rock layers for the compression and shear constituents, respectively. 

The real solid medium usually is thin-layered which at the definite assumption for the separate 
thick layer can be regarded as the gradient (or vertically-inhomogeneous) medium. Then the 

coefficient of attenuation [a'' " (m )] , characterizing the effects of energy dispersion in thin-layered 

beds, can be calculated statistically by means of the least squares method according to: 
a '• " = [ln(X + 2// ) j2 -ln(X + 2ii) j} ]/ AH j , where (A + 2/t ) } - the values of the longitudinal module of 

the elasticity on the upper ( / 7 ) and lower ( f 2 ) boundaries of the / - layer with the thickness ( AH j ) . 
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Simultaneously for receipt of bedded coefficients of the absorption \a'-"(ai )] , the interval 

values of the absorption coefficients, reduced earlier to bedded conditions, statistically average out by 
the least squares method on the depths intervals of layers respectively for the compression and shear 
constituents of the direct longitudinal wave. 

As the calculation of the e stimat e s of attenuation estimates (a'") and the absorption 

estimates [a' "(cJ)] is performed, using the independent data, that becomes possible owing to the 
numeral definition of the coefficient characterizing scattering of the longitudinal wave on the local 
intrastratal inhomogeneities according to the expression: cc ,,n disp -a f - n eff att -[a' " + a' "(aj )] , where 

5 r, " disp - the dispersion coefficient of disp e rsion on the local inhomogeneities for the j - layer with 

the thickness (AHj). The estimated coefficient of dispersion (5''" disp J are indicators of the 

stressed state of the separated layers and allow to obtain the additional information, for example, 
about the zones with the anomalous bedded pressures. 

The interpretating stage envisages the choice of the necessary totality of the measured oil- 
field parameters which with the sufficient completeness characterizes the reservoir and fluid- 
saturated properties of the deposit, moreover their information density guarantee the exact 
identification of the images of the oil-and-gas deposits. 

The totality of the significant oil-field parameters, including the coefficients of porosity ( <P) , 
the linear density of fractures (L d f ), granular (K g ) and fissured (K f ) permeabilities, residual water 

saturation (k' " rw ) and fluid saturation (T* " f s ) 9 quite unequivocally characterizes the capacitative, 

filtrational and fluid-saturated properties of the oil-and-gas-bearing layers formative the oil-and-gas 
deposit, moreover the above-mentioned parameters can be used as the indicators for the identification of 
the imag e of th e oil-and-gas deposit image . For the transformation of the petrophysical parameters to the 
form of the oil-field geology data - attributes of recognizable object the author offers the functional 
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relationship which are taking into account the termodvna mieal- thermodynamical conditions of the 
lay e rs of th e rocks layers . 

For the model of the pore-containing rock (or liquid-porous element of the deformable solid 
medium) granular permeability of the rock depends on the values of the average radius of the pore 
capillaries^, which can be calculated, as follows: R 0 =2{2D 0 £ /kT , where Q - the atomic 
volume; D 0 - the diffusion constant of atoms; k - Boltzmann's constant; T - Kelvin's temperature 
(°K). According to Pauzel's and Darsy's laws the granular permeability is defined following the 
formula: K g = 0 g R o 2 /8P 2 , where P - the meanderingness of pore channels coefficient and & g - the 

coefficient of granular porosity coefficien t. 

The granular porosity coefficient for the rock can be calculated as: 0 g =(p s -p)/(p s -p l ), 

where p,p s ,Pi are the common, solid and liquid phases of density, moreover the parameters p,p s ,Pi 

are defined, using the density log data, the lithological data and calibrated values of the absorption 
coefficients for the compression constituent of longitudinal wave. Hence, coefficient of the granular 
permeability, taking into account the previous formula, is as follows: K g = (0 g OD o X )/kTP 2 , where 

the parameter P is calculated by means of formula: P = 0~ m , moreover m = 1.82 (for the 

terrigenous rocks) and m = 2.03 (for the carbonate rocks). Then the residual water saturation 
coefficient for the porous medium can be defined according to: k' rw = (0 1/2 t , )/(2K g ) I/2 (l-0 g )P 9 

where r' - the average thickness of the lamina of the fixed water. 

For the model of the fracture-containing rock (or the fracture-brittle element of the deformable 
solid medium) the subhorizontal fissuring caused by the microfissures of the limited length, is 
characterized by means of such parameter as the linear density of fractures ratio, which is calculated 
here according to: L d f -(l t 25kTl)l DDt] , where / - the average length of microfissures; D- the 

diffusion constant of vacancies. Then the fissured porosity coefficient is calculated, using the formula: 
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0 f =L df b f , where b f - the average width (or opening) of microfissures. The opening fractures is 
calculated according to: b f = b Q exp[p f (p bed - p min )\, where b 0 - the initial width of microfractures; b f .- 
the co e ffici e nt of fissured compressibility coefficient ; p bedf p min> - the bedded and lateral mining pressures, 
accordingly respectively . At the same time the fissured permeability coefficient can be represented as 
follows: K f =(L d f b/*W~ 6 )/ 12 . Hence, the residual water saturation coefficient is defined according 

to: k* rw =(0 f 1/2 r" )/(3K f ) 1/2 (l-<P f ) y where z" - the average thickness of the lamina of the fixed 
water. 

Finally, the fluid saturation coefficients of the rock are calculated according to the formula: 
k f '" f s =(l-k' " rw ). If for the obtaining of the parametrical estimates of the reservoir properties 

(such as the permeabilities coefficients K g and K f ), the rheological parameters values (such as the 

coefficients of dynamical viscosity A* and rj ) are calculated previously, then the values of the 
parameters T 9 D o9 D 9 b f9 p bedm ,p min . > r' must be obtained additionally. So, for example, Kelvin's 
temperatures (T ) defines out of the temperature log data, and for the calculation of the opening of the 
microfissures the oil-field data concerning bedded ( p bed ) and mining ( p min ) pressures are used. 

The totality of the above-mentioned parametrical estimates for the determining of the imag e of th e 
oil-and-gas deposit image guarantees the correct identification of the images of various oil-and-gas 
deposits. After the r e c e iption reception of the estimates of the oil-field parameters procedures of the 
identification of the oil-and-gas deposit isimpl e m e nt e d is implemented b y means of the optimal deciding 
rules based on the Neyman - Pearson's criterion. 

The example of the experimental tests and the applied use of the declared method. For the 
borehole seismic investigations such technical devices, as digital recorders, explosive sources, 
three-component one- modular l evel sondes and other equipment have been used. The methods of 
the seismic data acquisition included: the excitation of the direct longitudinal wave by the near- 
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surface source of the oscillations (the one shotpoint), digital registration of the longitudinal wave 
using downhol e borehole and check three-component one modular leve l sondes and detailed 
observations (AH = 2.5-5 m) of the 3C VSP for the recording the displacement vectors for 
longitudinal and shear (or transverse) waves on the perspective in the productive relation depth 
interval in the investigated borehole. 

The results of the experimental diagnosis of the dynamical, rheological, absorbing, reservoir 
and fluid- properties of the rocks according to the detailed seismic investigation data, received on 
perspective in the productive relation an interval of depths {5030-5130} m in exploratory borehole of 
Lopushna - 13, which is located in the underthrust zone of the Pokutsko-Bukovynski Carpathians 
(Ukraine), are shown on Fig. (1 - 6). 

The geological problem which has been formulated before the seismic survey in the 
borehole Lopushna- 13 consisted in obtaining the quantitative estimates of the reservoir and fluid- 
saturated parameters for the possible perspective in the oil-and-gas-bearing relation thickness of 
Cretaceous depositions (argillites), which is located at the depths £5092.5 - 5105) m. The cap rock is 
located above the latter (low-porosity and low-permeable limestone). The results of the borehole 
seismic survey for the argillaceous thickness have shown the relatively high values of: granular 
porosity - 0 g = (3.4 - 8.8) %; Assuring - L d f = (48 - 65) m" 1 ; granular permeability - K g = (49 - 

202)*10" 3 mcm 2 and fissured permeability - K f = (16 - 40)*10" 3 mcm 2 , which was accompanied by 

of the most part of the layer with such values of the fluid saturation coefficients: k''" f s = 

(91-93)% (Fig.5-6). 

Thereby in the upper part of the layer (at depth of 5092.5 m) the zone of the tectonic 
dislocation has been diagnosed which is distinguished by the values of the interval coefficients of 
absorption a' "(c5) and the quality factors Q' "((d). Evidently the tectonic dislocation has changed 
sealing the deposit, earlier existing here, as the result of which only the residual oil in the upper part 
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of the layer [at the depths of (5092.5 - 5095)m] has been preserved. The subsequent oil - field tests 
of the given interval of depths have proved the above - mentioned prediction. 

Owing to the detailed investigations by the seismic survey method in the borehole of 
Lopushna - 13 the other zones with the anomalous values of the rheological and absorbing parameters 
of the rocks also have been identified (Fig. 3 - 4). According to the totality of the significant oil-field 
parameters (necessary for the identification of the imag e of th e oil-and-gas deposit imag e) at the 
depths of (5030 - 5065) m perspective layers have been detected which consist from the terrigenous 
rocks of the Neogene. 

Here the inhomogeneous reservoir is presented by the following parameters: 
0 g = (1.3 - 5.7) %; L d f = (19 - 44) m" 1 ; K g ={13 - 92.6)* 10" 3 mem 2 ; K f = (3.5 - 5 7. 9)* 10' 3 mem 2 

and k' r w = (6 - 16) % (Fig.5 - 6). The oil - and - gas deposit has the trilaminar disposition of fluid, 
which fills the pores and fractures of the rocks. In particular, the reservoir at the interval depths of: 
{5030 - 5037.5) m is gas-bearing; £5037.5 - 5047.5) m is oil-bearing; {5047.5 - 5052.5) m is water- 
bearing when the coefficients of fluid saturation vary within the following limits: from 84% up to 
94% (for the granular part) and from 92% up to 95% (for the fissured part). It is necessary to 
indicate, that the above - mentioned oil deposit has been missed earlier when only data of the 
standard complex of the well logging methods were used. 

The m e asur e d bedded coefficients of the effective attenuation (S 0 9,f ^ att ) 9 absorption 

\a'-"(a>) )* attenuation fa' " ) and dispersion (a t n disp ) hav e th e n e gativ e numerical valu e s in th e 

depth int e rvals: 5030 5050 m, 5092,5 5105 m and 5117,5 5030 m. The latest characterizes the pr e s e nc e 
of the zones with th e anomalously high pressur e in th e abov e mention e d lay e rs which is prov e d by the 
direct oil fi e ld for the compression and shear constituents of longitudinal wave, which are compared 
with lithology and data of the field measurements for bedded pressures and temperature s. Th e num e rical 
e stimates of th e bedd e d co e fficients of att e nuation, absorption and disp e rsion, r e c e iv e d for th e 
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r e sp e ctiv e compon e nt of th e direct longitudinal wav e , (the borehole of Lopushna -13), are represented 
on the Table 1 . 
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Table 1. 



Depth 
H[m] 


The compression constituent 
*10" 3 [m 1 ] 


The shear constituent 
*10° [m* ! ] 


Lithology 


The bedded 
pressure & 
temperature 


& cfl'.atl. 


a' (co) 


a ' 


G 'disp. 

r 


CC "cJT.att. 


a"(G> ) 


a » 


& "disp. 


5030- 
5050 


- 18,2 


6,7 
oil 


14,1 


-39,0 


-24,3 


0,9 


14,1 


-39,3 


Aleurolite 


7,8 MPa 
106°C 


5055- 
5065 


33,4 


2,9 
water 


17,4 


13,1 


32,1 


2,4 


17,4 


11,9 


Sandstone, 
argil lite 




5070- 
5090 


18,3 


2,8 
water 


1,2 


14,3 


11,5 


1,3 


1,2 


9,0 


Limestone 




5092,5- 
5105,0 


-38,0 


4,5 
water,oiI 


2,3 


-44,8 


-40,8 


1,6 


2,3 


-44,7 


Argillite 


7,6 MPa 
106°C 


5107,5- 
5115,0 


32,5 


9,7 
gas 


-3,2 


26,0 


23,5 


0,7 


-3,2 


26,0 


Clayey 
sandstone 




5117,5- 
5130,0 


-44,1 


2,3 
water 


2,4 


-48,8 


-44,2 


1,5 


2,4 


-48,1 


Dense 
limestone 


7,8 MPa 
108°C 



The measured coefficients of the effective attenuation ( 3''"^^ ) and dispersion 



( a'' " di ) have the negative numerical values in the depth intervals of: (5030-5050) rru (5092,5- 



5105) m and (5117,5-5030) m. The latest characterizes the presence of the zones with the 
anomalously high pressure in the above-mentioned layers, which is proved by the direct oil-field 
measurements of the bedded pressures and temperatures. 

Then, using the totality of the reservoir and fluid-saturated parameters of the rocks, the 
interval of depths interval of ( 5035-5050) m in the borehole of Lopushna-13 was recommended for 
the baring of the-productive layers. The tests have passed successfully with the obtaining of the oil 
influx with the initial flow rate of 0.8 m 3 /day out of the borehole that before had been considered as 
nonproductive. 

The results of the experimental tests allow us to predict the areas of the applied utilization of 
the declared method of wave diagnostics of the oil-and-gas deposit, in particular: 

- during the exploration drilling when it is necessary to detect the reservoir, cap rocks and 
the presence of the oil - and - gas - bearing layers; 
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- in the development drilling during the reservoir engineering of the oil-and-gas deposits 
when it is necessary to define more precisely the reservoir parameters of the productive layers, the 
quality of cap rocks, quantitatively to estimate the hydrocarbon saturation and to determine the 
locations of the gas - oil - water contacts; 

for the monitoring of the reservoir engineering of the oil-and-gas deposits. 
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CLAIMS 
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ABSTRACT 

The A method of wave diagnostics of the oil-and-gas deposit belongs to tho ar e a of the for 
seismic survey of th e oil and gas deposits on the basis of the complex of waves of the different 
types. It is performed by means of the insonification of the vicinity of the borehole by the seismic 
waves, the recording of the displacement vectors of the direct longitudinal and shear waves using 
the three-component seismic profiling in the investigated borehole, the data processing of the 
detailed profiling for the obtaining of the dynamical, rheological and absorbing parameters of the 
rock, and the geological interpretation of the above-mentioned parameters for the obtaining of the 
quantitative estimates of the reservoir and fluid-saturated properties of the rocks for each 
observation point lengthways depth of the borehole. Th e accuracy of the e s timat e d values of the 
dynamical param e t e rs which ar e r e c e iv e d by m e ans of th e param e trical analysis and d e coding of th e 
monotyp e t e mporal signals for th e compr e ssion and sh e ar components of th e dir e ct longitudinal 
wav e is ch e cked, using th e comput e r modelling of s e ismic signals. By m e ans of th e coproc e ssing of 
th e observation downhole and check data dynamical param e t e rs for th e downhol e data ar e 
corr e cted, moreover the changes of th e conditions of the wav e excitation and th e pr e s e nc e of th e 
filtrational e ffect from th e cov e ring thickness of th e rocks for e ach obs e rvation point lengthways 
d e pth of th e bor e hol e ar e taking into account that allows to calculat e th e dynamical param e t e rs of 
the impuls e respons e s and then to find th e r e sp e ctiv e p e trophysical param e t e rs of th e layers of the 
rocks. Furth e r th e p e trophysical param e t e rs ar e conv e rt e d into th e shape of th e oil fi e ld data th e 
coefficients of porosity, th e linear density of fractures, granular and fi s sured p e rm e ability, r e sidual 
wat e r saturation and fluid saturation of th e rook lay e rs, and calibrated valu e s of th e absorption 
coeffici e nts for the compr e ssion compon e nt of th e dir e ct longitudinal wav e ar e us e d as th e 
indicators of th e cont e nt of th e fluid and its typ e (oil, gas or wat e r). Diagnosing th e oil and gas 
deposit is r e alized by means of the us e of the totality of th e abov e m e ntioned oil fi e ld param e t e rs, 
moreov e r their information d e nsity guarantees the exact id e ntification of th e image of th e oil and 



37 



Marked Copy of Substitute Specification 
Application No. 10/649,941 
Attorney Docket No. 25585 

gas d e posit, and at tho sam e tim e it guarant ee s tho roceiption of th e r e liabl e information about th e 
e ff e ctiv e thickn e ss of layers and th e location of the gas oil wat e r contacts which can b e us e d for th e 
calculations of th e oil and gas r e serv e s. Th e pr e dicted ar e as of the appli e d utilization of th e d e clar e d 
m e thod of wav e diagnostics of th e oil and gas deposit: during th e e xploration drilling wh e n it is 
n e c e ssary to d e t e ct tho res e rvoir, th e cap rocks and tho pr e senc e of th e oil and gas b e aring lay e rs; in 
th e d e v e lopm e nt drilling during th e res e rvoir e ngine e ring of the oil and gas d e posits; for th e 
monitoring of th e res e rvoir e ngin ee ring of th e oil and gas deposits. 
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